Relation between inelastic electron tunneling and vibrational excitation of single 

adsorbates on metal surfaces 
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We analyse theoretically a relation between the vibrational generation rate of a single adsorbate by 
tunneling electrons and the inelastic tunneling (IET) current in scanning tunneling microscope, and 
the influence of the vibrational excitations on the rate of adsorbate motions. Special attention is paid 
to the effects of finite lifetime of the vibrational excitations. We show that in the vicinity and below 
the IET threshold the rate of adsorbate motion deviates from a simple power-law dependence on the 
bias voltage due to the effects of bath temperature and adsorbate vibrational lifetime broadenings. 
The temperature broadening appears to be confined near the threshold voltage within a narrow 
region of several fcsT, whereas the lifetime broadening manifests itself in a much wider region of 
applied voltages below the IET threshold. 
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PACS numbers: 68.37.EF, 68.43.Pq 



I. INTRODUCTION 



Since a pioneering demonstration of a single atomic 
switch by Eigler et al.fi in which high and low current 
state are realized via controlling Xe atom transfer be- 
tween a Ni substrate and a tip of scanning tunneling mi- 
croscope (STM), there have been very exciting develop- 
ments in the unique ability of STM to induce the motions 
and reactions of single adsorbates at surfaces. Exam- 
ples of such novel experiments are dissociation of decabo- 
rane on Si(lll)^ desorption of hydrogen from hydrogen- 
terminated Si(lOO)^ a step-by-step control of chemical 
reaction to form a biphenyl molecule from a iodoben- 
zene on Cu(lll)^ rotation and dissociation of an oxygen 
molecule on Pt(lll)^ rotation of acetylene molecule on 
Cu (100), 6 hopping of CO on Pd(llO), 7 and hopping and 
desorption of ammonia on Cu(100)i& 

The theoretical understanding of the physical mecha- 
nisms behind these motions and reactions of single ad- 
sorbates has made a slow but steady progress toward a 
full understanding of how tunneling electrons couple to 
nuclear motion of an adsorbate to overcome the poten- 
tial barrier along the relevant reaction coordinate. The 
power-law dependence of the atom transfer rate as a func- 
tion of applied voltage or tunneling current observed for 
the Eigler switch has been modelled as a potential bar- 
rier crossing between the potential wells formed by the 
interaction of the the adatom with the tip and the sub- 
strate, respectively. The atom overcomes the potential 
barrier through the stepwise vibrational ladder climbing 
by inelastic tunneling electrons. The so-called vibrational 
heating due to inelastic tunneling electrons has been pro- 
posed by Gao et al.fi Walkup et al.^ and Brandbyge 
& HedegarcUi [see also a detailed discussion in Ref. [l2j . 
Arrhenius-like expression for the rate of motion is charac- 
terized by an effective temperature T e g in the presence of 
inelastic tunneling current, otherwise a vibrational mode 
is in thermal equilibrium with a substrate temperature 
T. 



Another possibility of coherent multistep vibrational 
excitation by a single electron has been proposed by 
Salam et a/., 13 in analogy to a mechanism of desorption 
induced by electronic transition developed by Gadzuk. 14 
The coherent multiple excitation mechanism was shown 
to be important at low tunneling currents, where the av- 
erage time between successive electrons is longer than 
the vibrational lifetime. Dissociation rate of single O2 
on Pt(lll), where the vibrational relaxation rate of the 
0-0 stretch mode due to electron-hole pair excitations in 
the substrate is much larger than the tunneling current^ 
has been satisfactory described by this coherent multiple 
excitation process. 5 " However, one should bear in mind 
that since both coherent and incoherent multiple exci- 
tation mechanisms show a power-law dependence over a 
limited range of the bias voltage or the tunneling current, 
they can not be simply distinguished. 

Recently we have developediSiiLiS a theory of vibra- 
tional tunneling spectroscopy of adsorbates on metal sur- 
faces using the Keldysh-Green's function methodi2*2°i for 
an adsorbate-induced resonance model proposed by Pers- 
son and Baratoff. 2 i This has enabled us to elucidate 
the elementary processes of electron transport via adsor- 
bate state, and to obtain general formulae for the elas- 
tic and inelastic electron tunneling (IET), and the effec- 
tive vibrational distribution function within the second- 
order perturbation over the electron-phonon interaction 
of adsorbates. In comparison with the earlier papers^Siii 
where the Keldysh formalism was used as well, we take 
into account self-consistently the effects of the finite vi- 
bration lifetime on the IET spectra as well as on the 
vibration generation rate. 

In this work a relation between the IET current and 
vibrational generation rate is studied using the cou- 
pled Dyson and kinetic equations for causal and statis- 
tical phonon Keldysh-Green's functions. The stationary 
nonequilibrium phonon distribution function character- 
ized by the current-driven vibrational generation rate and 
the vibrational relaxation rate allows us to introduce the 
effective temperature and to describe self-consistently the 
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IET-induced vibrational heating. A primary interest is 
focused on the vicinity and below the threshold of the 
vibrational exc itation. This situation was not investi- 
gated in Refs. lijl lull lll2t where the limit of large applied 
voltage was considered, eV 3> CI (CI is the vibrational 
mode energy). As we show, the dependencies of the vi- 
brational generation rate and of the multi-step adsorbate 
excitation rates on the bias voltage deviate from simple 
power laws due to the effects of bath temperature and a 
finite lifetime of the vibrational mode excited by tunnel- 
ing electrons. Whereas the bath temperature broadening 
is restricted to a narrow region within several fcgT, the 
vibrational lifetime broadening causes a strong overheat- 
ing effect in a much wider bias voltage region below the 
IET threshold. 

The paper is organized as follows. In Section ITU the 
general formulation of the problem in terms of Keldysh- 
Green's functions is outlined, and the kinetic equation for 
vibrational excitation is derived. In Section IlIII we ana- 
lyze the expressions for the IET vibrational generation 
rates and IET current. The physical meaning of IET vi- 
brational generation rate is discussed in Subsection lHI Al 
In subsection llll Bl the limiting cases of low bath temper- 
ature and no vibrational broadening are analyzed, and 
it is shown that in this limit our results coincide with 
those received previously in Refs. l9l lTollTlllT^ . Numerical 
examples are discussed in Subsection llll CI 



II. THEORY 

We consider the same Hamiltonian as used be- 

£ ^12.16.17.18.21 

H = ekcjck + e v c\c v + e a c\c a + Cl(tfb + 1/2) 
fe p 

^ ] (Vpa,CpC a + h.C.) 
fe p 

+ x4c a (6t+&), (1) 

where the energies and annihilation operators of a sub- 
strate, a tip, an adsorbate orbital and a vibrational mode 
are denoted by e p , Ek, s a , CI and c p , Ck, c a , 6, respectively. 
The tip and substrate systems are assumed to be in 
thermal equilibrium at the same temperature T, and to 
have independent chemical potentials fit and fi s , respec- 
tively, whose difference corresponds to the bias voltage 
eV = n s — fit- Electronic tunneling matrix elements Vk a 
(tip-adsorbate) and V pa (substrate-adsorbate) give rise 
to stationary current between the tip and the substrate 
through the adsorbate orbital. The distribution function 
for electrons of the tip and substrate systems, n t and n s 
respectively, are given by the Fermi distribution function. 
X is the vibration- adsorbate orbital coupling constant. 

In the Keldysh-Green's function method 19 the 
nonequilibrium process of the electron tunneling be- 
tween the tip and substrate through the adsorbate level 



is described by the coupled Dyson and kinetic equa- 
tions for causal and statistical Keldysh-Green's func- 
tions. This coupled description of the dynamical and 
statistical properties of the nonequilibrium system al- 
lows, within the Keldysh-Green's function method, to 
find s el] -consistently the spectra and occupation num- 
bers changes which are introduced into the system by 
the external action (the tip-substrate bias voltage in our 
case), mediated by the interactions in the system. In 
what follows, the implementation of the the Keldysh- 
Green's function method for the Hamiltonian is out- 
lined by Eqs. (|2*|)- (|10f> without detailed explanations; the 
latter can be found in the original paper by Keldyshi 9 -, 
see also in Refs. 

The kinetic equation for vibrational excitations 
(phonons) takes the form 



dN, 



P h 



at 



= J g[lI+-(a;)D-+H-n-+(a;)D+-( W )]. 



(2) 

Here the phonon Keldysh-Green's functions D ±T are 
given by 



D + -(lo) = n ph (uj)[D r (u)-D a (uj)} 

= -2iirn ph (uj)p ph (Lj), 

D-+{uj) = [l + n ph H] 

= -2wr[l + n p h(w)]/9ph(w), 



(3) 
(4) 



where n p h(w) is the vibration occupation function to be 
determined self-consistently from the kinetic equation, 
and the phonon retarded and advanced Green's functions 
are given by 



1 



oj - 0+ i7 e /i(w)/2' 
1 

to - CI - ij eh (uj)/2' 



(5) 



Here 7 e / l (w) = — 2ImH r (where n r is the retarded 
phonon polarization operator) is the vibrational damping 
rate due to electron-hole pair excitations in the substrate 
and tip^ We omit here the real part of n r , which gives 
a red shift of the resonant frequency CI. 

The phonon polarization operators an d n r are 

calculated in the second order perturbation over x- F° r 
example, 



rr 



J ^G + a -{e + u)G- + (e) 



(6) 



n r =^ / ^[G^e)F a ( e + U )+F a {e)G r a (e + u;)], 

(7) 

where the retarded, advanced and statistical Keldysh- 
Green's functions of the adsorbate are connected as 

F a (e) = [l-2n a (e)}[G r a (s)-G a a (e)] (8) 
= -2i?r[l-2n a (e)}p a (e), 
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with the adsorbate density of states p a and occupation 
function n a , 



p a (e) = (l/n)A/[(e-e a ) 2 + A 2 }, 



n a (s) 



n s (e)A s + n t (e)A t 



(9) 



(10) 



partial width due to the hybridizations between the tip 
and the adsorbate (substrate and adsorbate). The cor- 
rections (due to the vibration-adsorbate interaction) to 
the adsorbate density of states and to the adsorbate occu- 
pation function, Eqs.© and (f 101) . respectively, are in the 
second order over the coupling constant \- Thus, they 
can be neglected in Eas. (|6l7f) which are already written 
in the second order over \- 



Here A ~ A t + A s is the resonant width of the adsor- 
bate orbital, A t{s) (e) = nJ2k( P ) \ v k( P )a\ 2 S(e - e k ) is the 



Then Eq. © leads to 



= 2lT X 2 I dedujp a (e)p a {e + uj)p ph (uj) [1 - n a {e)\ n a (s + w) - j duj-/ eh (uj)n ph (uj)pp h (uj), 



(11) 



where /o p h(e) = (l/27r)7 e /i(w)/[(ti> — £1) 2 + 7^(^)/4] is the vibrational density of states with the broadening, and the 
vibrational damping rate (or linewidth) j e h due to electron-hole pair excitation is 



leh{u) = 2irx 2 / de[n a (e) - n a {e + u)} p a (e)p a {s + u) w 2ttx 



2 A s P 2 a (Ps) + A t p 2 a {p t 
A 



(12) 



r 



The last approximate equation in (|12[) is valid in the limit 
of low-temperature, ksT <C lu and slow- varying p a (s) 
over eV. 

Adding and subtracting integral term 
/ duj^f e h(uj )nr(w)/3 p h(w) to the kinetic equation (| 1 1 p . it 



can be transferred to a more transparent form, showing 
that in the absence of the applied voltage the stationary 
phonon distribution function is the equilibrium one, 
n T {uj) = ( e "/ fc * T - l)-\ 



n = 9N ph 
dt 

Here 

.A, A* 



J dwpp^uj) {T in (uj) - 7 e h(w) [n-ph(^) - n T (oj)]} . (13) 



r in (w) = 2ir\ . 2 / dep a (e)p a (e + u) [n s (e + u) - n t (e + cj)\ [n s (e) - n t (e)] , (14) 



is the vibrational generation rate function. 



III. RATES OF VIBRATIONAL GENERATION 
AND ADSORBATE MOTION BY TUNNELING 
CURRENT 

Kinetic equation Eq. I|13|) in the stationary regime 
allows us to find the stationary distribution function 
n p h(w) which deviates from the equilibrium one nr(w). 
The consequences of this deviation on the rate of adsor- 
bate motion can be analysed within the truncated oscil- 
lator model, introduced and fully described in Ref. 0. 
The main quantities in the truncated oscillator model** 
are the the vibrational excitation and deexcitation rates 
between the nearest-neighboring levels of the harmonic 
oscillator, Tf and respectively. Knowing it be- 



comes possible, as shown in Ref. 0, to introduce the 
adsorbate effective temperature 

T -=£( ln ft)"= < 15 > 

and, solving the oscillator master equation, to estimate 
the rate of adsorbate motion 

fl(v) = „r t (Et)«-\ (i6) 
1 1 

where n is the number of the vibrational levels within 
the truncated harmonic approximation for the potential 
well, and (n — is assumed to be close to the barrier 
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height. Equation (|16fl can be also written^ as Arrhenius- 
like expression via the effective temperature, 



R(V) = nr T exp[- 



(n- i)n 
k B T eS (ny 



(17) 



In what follows, we generalize the equations for the 
excitation and deexcitation rates from Ref. I12I in order 
to include the effect of the finite lifetime broadening of 
the vibrational excitations. The solution of the kinetic 
equation 113|) for the stationary vibrational occupation 
number 



N, 



P h 



dwp p h(w)n p h(w) 



(18) 



can be read as the effective Bose-Einstein distribution 
function 



iVph = [exp(n/fc B T c ff) - 1]" 



(19) 



in the presence of the IET-driven generation T^t in com- 
peting to its damping 7. Comparing the latter with the 
definition of the effective temperature via the excitation 
and deexcitation rates, Eq. I|15fl . arid taking into account 
the explicit form of the kinetic equation (|13fl , these rates 
can be written as 

r T = n T (n)% h + r iet , r x = [n T (0) + \]% h + r ict , (20) 

where, in contrast to Ref.lT^. the effect of the vibrational 
lifetime broadening is included into the generalized form 
of the IET vibration generation rate via the integration of 
Tin over the lifetime-broadened vibrational density func- 
tion, 



and 



7eh 



dw/9 p h(w)r in (w), 



/ dm p ph (uj)"f eh (u)n T (lj) 



(21) 



(22) 



is an averaged vibrational linewidth. Actually, because 
in reality the vibrational linewidth is always relatively 
small, 7 e /j <C and 7 e /i(w) is a slowly changing function 
at u> ~ f2, a good approximation for Eq. (|22|l is j e h ~ 

7eh.(^)- 

We would like to note here that, in principle, we have 
to include additively (see the discussion above) the relax- 
ation rates via other channels due to phonon excitations 
in the substrate 7 P h, if they are important, and j e h is 
replaced by total relaxation rate, 7 = jeh + 7 P h^ 



A. Physical meaning of Fi ct 

Let us discuss the physical meaning of rj ot from the 
point of view of the elementary processes of the elastic 
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FIG. 1: Schematic representation of the inelastic electron 
transfers through the adsorbate from tip to substrate with 
emission and absorption of phonons Rf 3 and R^~ s , (a) and (b), 
respectively. Same from substrate to tip with emission and 
absorption of phonons R+ t and Rj t , (c) and (d), respectively. 



and inelastic electron tunneling through the adsorbate- 
induced resonance. 

As has been noted in Ref. lljj, the main difference be- 
tween the elementary processes composing the inelastic 
current and the vibrational generation rate, is the fact 
that at T ^ there are inverse (hole) inelastic currents, 
which come with different signs to the total current and 
to the vibrational generation rate. The expression for the 
inelastic tunneling rate in the second order over Xi [ see 
Eqs. (14) from Ref. 17.], can be written as 



.A, A* 



(23) 



- 27TX 2 



A 2 



de / dcop^(e)p^(e + Lu)p ph (Lu) 



x {-n s (e + Ld)[l-n t (e)}[l + n ph (Lj)} 

- n s (e) [1 -n t (e + u)]n ph (uj) 

+ [l-n s (e)]n t (e + uj)[l + n ph (uj)} 

+ [1 - n s (s + u)] n t (e)n ph (uj) } . 

As schematically illustrated in Fig. 1 each term of this 
equation has its own inelastic tunneling process, and it 
can be deconvoluted into 



Rin - +R st + R s 



Rts 



r; 



(24) 



Here R st denotes the electron transfer rate through 
the adsorbate level with vibrational emission/ absorption 
from substrate to tip (or a hole from tip to substrate), 
and Rf s is the reverse process from tip to substrate. In 
Figure 1, the initial electron states are shown as empty 
circles below the corresponding Fermi levels. The final 
states are shown as solid circles. In order that the par- 
tial currents Rf^ are nonzero, the final states should be 
empty. For example, at T — and the applied voltage 
shown in Fig. 1, the currents shown in panels (c) and 
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(d) are zero. They become slightly nonzero at finite tem- 
peratures, due to thermal excitation of holes below the 
Fermi levels. 

We turn now to the calculation of the vibrational gen- 
eration rate. On the one side, it is zero in the station- 
ary situation, Eq. (|13|1 . On the other side, it can be 
calculated within the same accuracy over \ as Eq (|24ll . 
summing up the electron transfer rates with vibrational 
generation and subtracting the rates with absorption, in- 
cluding the processes of generation and absorption leav- 
ing the electrons in substrate and in tip. The latter do not 
contribute into current and, thus, are absent in Eq (|24() . 
Then, 
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where, e.g., 



TDSt _ + ) _ p(-) , p( + ) _ p(-) 



(26) 



compare the signs here with that in Eq. I|24(l . After some 
algebra we obtain, 



<h = < h {n ph (w)} 



(27) 



= 2**^ y y ^p(°)( £ )pi )( e + c)p ph ( w ) 

x (n s (e + w) [1 -n t {e)] + n t (e + to) [1 -n a (e)] (28) 
+ n ph (w) [n«(e + lo) + n t (e + uS) - n s (e) - n t (e)]). 

It appears that, in order to calculate r; Gt , Eq.JSlJ, via 
Rp h , we have to replace in Eq. I|28|l the stationary 
nonequilibrium distribution function n p h(w) by the equi- 
librium one. Using notation of Eq. 1|27|) . it can be written 

as 



r,, 



(29) 
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FIG. 2: Voltage dependences of the inelastic current 7?i n and 
vibrational generation rate Ti et (lines and lines with circles, 
respectively) , calculated for T — 4K and different vibrational 
lifetimes (a); for r e h = 2 ps and different temperatures (b). 
Fhe vibrational frequency is taken Q = 100 meV. See other 
parameters in the text. 



This means that there is a small but systematic differ- 
ence between r; et Eqs. I|29l2ti[) and absolute value of IET 
rate |i?; n | Eqs. Q23I24J1 . growing with the increase of de- 
viation of stationary vibration population n p h from the 
equilibrium one tit, see also Eq. (1321) below. 



B. Limit of T = and no vibrational broadening 

Before presenting the results of numerical calculation 
of ri C t,-Rin, and R at finite temperatures and including 
the effect of vibrational broadening, it is helpful to ana- 
lyze the situation in the limit \eV\ ^> f2 ksT, neglect- 
ing the vibrational broadening, p p h (u) = 8{u - fl). For 
simplicity, we will also assume A s ^> At. We expect to 
obtain the results equivalent to those in Gao et alH 

Equation (|18|l then reduces to 



r,, 



lenity ' 



(30) 



with 

Tiet 



2 A s A t F(x) 



A 2 



Ai 7eh(n) ~cr- 



(31) 



where F(x) — xOix), x — \eV\ — Q, and rj — sign(F) is 
the sign function and 9 is the step function. In the last 
step we have used p a ipt — vty — Pai.Pt)- In the limit 
\eV\ 3> Vi Eq. 1|31[) becomes a linear function of applied 
voltage, in agrement with the result obtained by Gao 



et 



from a Boltzmann distribution among different 



vibrational levels. 

As to a relation between i?; n and T^t in this limit, 
comparing Eqs. (|23|l and (|29|1 . we find that 



\Rn 



leh- 



i\eV\ - Q)\eV\ 



Qi\eV\ 



n). 

(32) 
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FIG. 3: Voltage dependences of the effective temperature, 
calculated for bath temperature T = 4 K and different vibra- 
tional lifetimes (a); for r e u = 2 ps and different temperatures 
(b). Short dotted line in both panels shows the effective tem- 
perature in the simple approximation neglecting the phonon 
broadening effect, r e h = oo, and at T = 0, ri e t(V) cx \eV\ —Q. 



FIG. 4: Voltage dependences of the logarithm of the trans- 
fer rate, calculated for bath temperature T = 4 K and 
different vibrational lifetimes (a); for r e /, = 2 ps and dif- 
ferent temperatures (b). Short dotted line in both panels 
shows the transfer rate in the simple approximation neglect- 
ing the phonon broadening effect, r e h ~ oo, and at T = 0, 
R(V) oc (|eV| - Q,) n , n = 3. 



Because (A t /A s ) 2 <C 1, |i? in « T iet . The linear relation 
between Ti ct (or approximately i? in ) and the bias voltage 
holds for eV 3> f2, which is a case of a Xe atom transfer 
between a Ni substrate and a STM tip (0=4 meV and 
V =20-200 mV)±±£ 

As to the adsorbate motion rate, in this limit a linear 
dependence of T ict (V) on \eV\ — O, (Eq. ||31J) and a con- 
dition of 7 > rj c t(T^) yield the power-law dependence 



R(V) oc (\eV\ - Cl) n (33) 

for Eq. (|16|l . in agreement with the analysis in Ref. 0. 
The opposite case 7 Ti ct (V) implies infinite effective 
temperature which is not consistent with an assumption 
of activation- type-law for a reaction rate Eq. (|17|l . 



C. Numerical examples 

We are in a position to compare the relation between 
inelastic tunneling rate Ri n ,e and vibrational generation 
rate Ti ct as a function of V for several values of the vi- 
brational relaxation time and temperature. Numerical 
calculations are performed for a set of parameters: £ Q =2 
eV above the Fermi level of the substrate, A s = 1.0 eV, 
A t =0.025 eV, ft =100 meV, and different temperatures 
T = 0.4, 1, 4, and 40 K. The electron- vibrational coupling 
constants were taken \ = 500,354,250 and 177 meV, 
which at eV = and given A Sj t correspond to vibra- 
tional linewidth 7^ = 0.66, 0.33, 0.16, and 0.08 meV 
(or vibrational lifetimes r e h — fr/jeh =1 ; 2, 4, and 8 ps), 
respectively. 

Figure [2] shows the calculated vibrational generation 
rates Eq. I|21|l (circles) and the corresponding inelastic 
currents Ea. (|23|l (lines) at fixed temperature for different 
vibrational lifetimes (panel a) and for fixed vibrational 
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lifetime for different temperatures (panel b). The gen- 
eration rates are given in nA and also in /xeV (left and 
right vertical axes in Fig. |5J). Because the vibrational 
linewidths are taken in the range of 0.088 - 0.66 meV, 
it is clear that all the calculated rates are well within 
the condition 7 3> Fi e t(V). It is seen that, in agreement 
with the analysis above, the inelastic current starts to 
exceed slightly the vibrational generation rate above the 
threshold of inelastic tunneling. As to the latter, it grows 
approximately linearly, cx \eV\ — SI above the threshold. 
The differences from the simple linear dependence arise 
near the threshold. Instead of abrupt switching on of the 
inelastic current and vibrational generation rate, there is 
a gradual growing of both. The physical reason is the 
temperature and lifetime broadening of the IET thresh- 
old. The temperature broadening has been analyzed pre- 
viously (see, e.g., in Ref l23T) . The lifetime broadening was 
not analyzed, at least in part of its influence on the over- 
heating and adsorbate motions. As we will see below, 
a different influence of temperature and lifetime broad- 
enings allows, in principle, to separate these two mech- 
anisms from the voltage dependence of the overheating 
and of motion rate below the IET threshold. 

It becomes more evident from Fig. where the calcu- 
lated voltage dependencies of the effective temperature 
Eq. (|15|l for the same parameters as in Fig. [21 are dis- 
played. 

Note that above the IET threshold, \eV\ > SI the over- 
heating does not depend on the vibrational lifetime as 
well as on the temperature and approaches a universal 
linear dependence Eq. (|3I|1 oc \eV\ — SI. The reason for 
this is the cancellation of % 2 in the nominator and de- 
nominator of Ti C t/j e h, if the e-h scattering Eq. (|12|) is 
the main channel of vibrational relaxation, 7 « "f e h . As 
a result, the stationary vibrational occupation number 
Q30JI. effective temperature ()I5|I. and rate of adsorbate 
motion l|16|l do not depend on the adsorbate-vibration 
coupling constant. We have already mentioned this fact 
previously^ This effect takes place only above the IET 
threshold. 

Whereas near and below the threshold, \eV\ < SI, the 
effective temperature becomes lifetime and bath temper- 
ature dependent. First of all note that, in a seeming 
contradiction with a simple understanding that a shorter 
lifetime means better cooling, an opposite dependence 
is seen in Fig. the overheating grows with the life- 
time decrease! The explanation is very simple: at low 
temperatures and below the IET threshold, the only rea- 
son of the nonzero inelastic current and overheating is 
the finite lifetime broadening. Which grows with the 
lifetime decrease. Mathematically, an additional power 
of x 2 is brought into r ict below the IET threshold be- 
cause of integration of r in (cj) oc \ 2 m Eq. (|2"T|l over the 
tail of the vibrational density of states at large detuning 
7ph(w) oc j e h/(uj — SI) 2 , which is oc \ 2 too. 

Figure |2h> shows that, as anticipated, the tempera- 
ture broadening is negligibly small in comparison with 
lifetime broadening, for lower temperatures fc^T < 7^. 



For example, all T e g curves in Fig. |2h>, calculated with 
7e/i = 0.33 meV, merge into a single limiting curve start- 
ing from bath temperatures T < 1 K (ksT < 0.9 meV). 
This limiting curve corresponds to a T = solution for 
effective temperature with a finite vibrational lifetime. 

Another interesting thing can be seen in Fig. the 
temperature broadening at higher temperatures, fc^T > 
7c/!, is localised within a finite range of several ksT near 
the IET threshold, whereas the finite lifetime broaden- 
ing is characterized by a long tail of strong overheating 
even at \eV\ <C ft. The explanation is that the latter is 
due to Lorenzian tail of the vibrational density of states, 
whereas the temperature broadening has a Gaussian tail. 

Figure 0] plots R{V) Ea. ((T7|l for the same set of pa- 
rameters as before and barrier hight equal to n ~ 3 vi- 
brational excitation quanta. It shows that although with 
the increase of voltage above the IET threshold eV > SI, 
the rates of motion reach eventually the limiting power 
law dependence Eq. (|33|l , in the vicinity of the threshold 
this dependence is different, depending on the relative 
values of fc^T and 7^. And below the threshold, at 
\eV\ <C SI — fcflT, the rate of motion is fully due to the 
vibrational lifetime broadening, thus making it to grow 
with the decrease of the vibrational lifetime. 



IV. CONCLUDING REMARKS 

We have studied a formal relation between the inelas- 
tic tunneling current, vibrational generation rate and vi- 
brational heating caused by the current-induced vibra- 
tional excitation of single adsorbates on metal surfaces. 
The nonequilibrium process of the electron tunneling be- 
tween the tip and substrate through the adsorbate level 
is described by coupled Dyson and kinetic equations of 
the Keldysh-Green's functions. The inelastic tunneling 
current and the vibrational generation rate are obtained 
within the second-order perturbation with respect to 
the electron-vibration coupling based on the adsorbate- 
induced resonance model. We show that the vibrational 
generation rate slightly differs from the absolute value of 
the inelastic tunneling rate. The linear dependence of the 
vibrational generation rate on the bias voltage is valid for 
the extreme cases of \eV\ 3> SI and the slowly varying ad- 
sorbate density of states over the applied voltage. Under 
these conditions the rate of adsorbate motions via the 
step- wise vibrational ladder climbing due to the inelastic 
tunneling current is characterized by a simple power-law 
dependence on the bias voltage or the tunneling current. 
It is found that the effects of bath temperature and vibra- 
tional damping manifest themselves in the vicinity and 
below the IET threshold, at \eV\ < SI. The temperature 
broadening manifests itself only at relatively large bath 
temperatures, ksT > 7^. The important thing is that 
whereas the bath temperature broadening is confined to 
a narrow range of several fcgT, the vibrational lifetime 
broadening effect extends to a much larger interval of 
voltages below the IET threshold. 



8 



ACKNOWLEDMENTS We would like to thank T. 
Mii for useful discussions. This work was financially sup- 
ported in part by the Russian Foundation for Basic Re- 
search and the Russian Ministry of Science. S.T. thanks 



Japan Society for the Promotion of Science fellowship 
program (JSPS) (No. S02227) for supporting his visit to 
Toyama University. H. U. is supported by Grant-in-Aid 
for Basic Research (B) (No. 15310071) from JSPS. 



1 D.M. Eigler, CP. Lutz and W.E. Rudge, Nature 352, 600 
(1991). 

2 G. Dujardin, R.E. Walkup and Ph. Avouris, Science 255, 
1232 (1992). 

3 T.C. Shen, C. Wang, G. Abeln, J.R. Tucker, J.W. Lyding, 
Ph. Avouris, R. Walkup, Science 268, 1590 (1995). 

4 S.-W. Hla, L. Bartels, G. Meyer and K.-H. Rieder, Phys. 
Rev. Lett. 85, 2777 (2000). 

5 B.C. Stipe, M.A. Rezaei, W. Ho, S. Gao, M. Persson and 
B.I. Lundqvist, Phys. Rev. Lett. 78, 4410 (1997). 

6 B.C. Stipe, M.A. Rezaei and W. Ho, Phys. Rev. Lett. 81, 
1263 (1998). 

7 T. Komeda, Y. Kim, M. Kawai, B.N.J. Persson and H. 
Ueba, Science 295, 2055 (2002); B.N.J. Persson and H. 
Ueba, Surf. Sci. 502-503, 18 (2002). 

8 J.I. Pascual, M. Lorente, Z. Song, H. Conrad and H.- 
P.Rust, Nature 423, 525 (2003). 

9 S. Gao, M. Persson and B.I. Lundqvist, Solid State Com- 
munications 84, 271 (1992); J. Electron Spectro. and Re- 
lated Pheno. 64/65, 665 (1993). 

10 R. E. Walkup, D. M. Newns and Ph. Avouris, J. Electron 
Spectro. and Related Pheno. 64/65, 523 (1993). 

11 M. Brandbyge and P. Hedegard, Phys. Rev. Lett. 72, 2919 
(1994). 



S. Gao, M. Persson and B.I. Lundqvist, Phys. Rev. B 55, 
4825 (1997). 

13 G.P. Salam, M. Persson and R.E Palmer, Phys. Rev. B49, 
10655 (1994). 

14 J.W. Gadzuk, Phys. Rev. B 44, 13466 (1991). 

15 B.N.J. Persson, Chem. Phys. Lett. 139, 457 (1987). 

16 S. Tikhodeev, N. Mingo, K. Makoshi, T. Mii and H. Ueba, 
Surf. Sci. 493, 63 (2001). 

17 T. Mii, S.G. Tikhodeev, and H. Ueba, Surf. Sci. 502-503, 
26 (2002). 

18 T. Mii, S.G. Tikhodeev, and H. Ueba, Phys. Rev. B 68, 
205406 (2003). 

19 L.V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964) [Sov. 
Phys. -JETP 20, 1018 (1965)]. 

20 P.I. Arseyev and N.S. Maslova, Zh. Eksp. Teor. Fiz. 102, 
1056 (1992) [Sov.Phys.-JETP 75, 575 (1992).] 

21 B.N.J. Persson and A. Baratoff, Phys. Rev. Lett. 59, 339 
(1987). 

22 If the electron-hole pair scattering is not the only mecha- 
nism of the phonon relaxation, other mechanisms should 
be included into Tl r . 

23 S. Gao, Surf. Sci. 313, 448 (1994). 



